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ABSTRACT 


The  report  presents  the  results  of  a  numerical  Investigation  of  a 
novel  recurrence  relation  for  iteratively  solving  modified  integral 
equations  for  calculating  free-surface  potential  flow.  Numerical 
results  are  presented  for  the  simple  case  of  wave  radiation  in  the 
zero-  and  infinite-frequency  limits. 

ADMINISTRATIVE  INFORMATION 

This  research  was  funded  under  the  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center's  Independent  Research  Program.  Program  Element  61152N,  Task 
Area  ZR0230101,  and  Work  Unit  1843-045. 

INTRODUCTION 

Most  of  the  numerical  methods  reported  in  the  ship-hydrodynamics  literature 
for  calculating  potential  flow  about  ship-hull  forms,  whether  free-surface 
effects  are  taken  into  account  or  are  neglected  as  in  the  so-called  double-hull 
approximation  (in  which  the  free  surface  is  treated  as  a  rigid  wall),  closely 
follow  the  pioneering  method  of  Hess  and  Smlth^.  Briefly,  this  well-known 
numerical  method  consists  in  determining  the  density  of  an  auxiliary  distribu¬ 
tion  of  sources  on  the  body  surface,  in  the  classical  manner  described  by 
Kellogg^,  by  numerically  solving  an  Integral  equation.  This  Integral  equation 
is  solved  by  approximating  the  surface  of  the  body  by  an  ensemble  of  planar 
panels  and  assuming  a  constant  (uniform)  source  density  within  each  panel.  The 
integral  equation  is  then  approximated  by  a  system  of  linear  equations,  which  is 
solved  by  Inverting  the  coefficient-matrix,  so-called  "matrix  of  influence 
coefficients". 

A  number  of  variations  and  extensions  of  the  method  of  Hess  and  Smith  have 
been  proposed.  In  particular,  alternative  integral  equations  have  been 


formulated  —  by  using  an  auxiliary  distribution  of  normal  dipoles  (instead  of 
sources  as  in  the  method  of  Hess  and  Smith)  or  by  directly  applying  Green's 
third  identity  to  the  velocity  potential  and  the  fundamental  source  potential 
(Green  function)  — *  and  solved  numerically.  Important  improvements  to  the  basic 
Hess  and  Smith  method  have  recently  been  incorporated  into  the  so-called  higher- 
order  panel  method^,  in  which  curved  panels  supporting  linear-source  and 
quadratic-dipole  distributions  are  used. 

An  alternative  approach  for  calculating  potential  flow  about  a  body  in  an 
unbounded  fluid,  in  tihlch  a  modified  integral  equation  for  determining  the 
velocity  potential  is  solved  iteratively,  has  been  proposed^ This  alter¬ 
native  approach  has  also  been  extended  to  the  calculation  of  free-surface  poten¬ 
tial  flow8^>^>^^.  A  main  recommendation  of  this  alternative  approach  is  that  the 
modified  Integral  equations  upon  which  the  approach  is  based  have  regular 
kernels,  whereas  the  traditional  Integral  equations  used  in  the  previously- 
mentioned  methoc 5  have  singular  kernels.  Another  potential  advantage  of  the 
alternative  approach  resides  in  the  use  of  an  iterative  solution  procedure, 
which  is  well  suited  to  the  analysis  of  complex  systems  requiring  a  large  number 
of  panels^.  This  report  presents  the  results  of  a  limited  numerical  Investiga¬ 
tion  of  the  novel  recurrence  relation  proposed  in  Noblesse^®.  More  precisely, 
this  recurrence  relation  is  studied  here  for  the  simple  case  of  wave  radiation 
in  the  zero-  and  infinite-frequency  limits. 

INTEGRAL  EQUATIONS 

This  study  considers  potential  flow  due  to  translation  or  rotation  of  a 
rigid,  nonlifting  body  through  an  incompressible  and  Invlscid  fluid.  The  sur¬ 
face  of  the  body  is  denoted  by  S,  and  n  represents  the  unit  normal  vector  to  S 
pointing  Inside  the  body.  Furthermore,  T  (x,  y,  z)  and  C  (5.  n,  C)  represent 
the  position  vectors,  with  respect  to  a  rectangular  system  of  coordinates 
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attached  to  the  noving  body,  of  any  two  points  on  the  body  surface  S.  The  coor- 
dlnates  x  and  $  are  adlaenslonal,  with  respect  to  some  length  characterising 
the  size  of  the  body  taken  as  reference.  The  velocity  potential,  ^  say,  of  the 
flow  likewise  Is  adlaenslonal,  with  the  body  characteristic  length  and  sone 
characteristic  speed  of  Its  motion  used  as  reference.  The  normal  component  of 
the  fluid  velocity  at  any  point  IT  on  the  body  surface,  given  by  -  V^*?,  Is 
equal  to  the  normal  component  of  the  velocity  of  S  at  x,  and  thus  Is  known  for 
all  six  fundamental  motions  of  the  body  (translation  and  rotation  about  the  x, 
y,  z  axes). 

The  Green  function  appropriate  for  an  unbounded  fluid  Is  given  by 
4xG(C,x)  ■  -1/r  %fhere  r  -  [(5-x)^  +  (n-y)^  +  Application  of  Green's 

third  Identity  to  the  velocity  potential  ^(x)  and  the  foregoing  free-space  Green 
function  G(^,x)  yields  the  following  classical  integral  equation  for  determining 
the  potential  at  any  point  C  on  S: 


♦  (0/2  -  f(0+L(c;^), 


where  f(^)  Is  the  known  (since  ♦{}  Is  specified  on  S)  potential  defined  as 


f(0  ■  ~/s  G(C,x)tn(x)da(x), 


and  L(^;^)  represents  the  Integral  transform  of  ♦  defined  as 


b(5;^)  ■  /s  Gn(5,x)^(x)da(x). 


In  equations  (2)  and  (3),  da(x)  represents  the  differential  element  of  area  of 
the  surface  S  at  the  point  IT. 

The  previously-mentioned  modified  Integral  equation  corresponding  to  the 
classical  integral  equation  (1)  takes  the  form^>^>^ 


♦  (5)  -  f(0+L'(5;^), 


where  f(S)  Is  the  known  source  potential  given  bjr  equation  (2)  and  L' is 
the  integral  transform  of  4  defined  as 

■  /s  ®n(T.x)[4(^-4^)lda(T).  (5) 


A  noteworthy  feature  of  the  modified  integral  equation  (4)  is  that  the  integrand 
Gn(5 ,x)[4(x)-4(5 )]  of  the  Integral  (5)  remains  finite  as  the  "Integration 
point"  X  approaches  the  "calculation  point"  C* 

The  radiation  potential  for  a  body  intersecting  a  free  surface  must  satisfy 
the  free-surface  boundary  condition  -  (w  L/g)  4  “  0  on  the  equilibrium 

plane  of  the  free  surface  taken  as  the  plane  z  •  0  with  the  z  axis  pointing 
upward,  %ihere  (o  >  radiant  frequency  of  oscillatory  motion  of  the  body,  L  *  body 
characteristic  length  and  g  >  acceleration  of  gravity.  In  the  zero-frequency 
limit  u  L/g  ■  0  and  the  infinite-frequency  limit  u  L/g  ■  •,  the  free-surface 
boundary  condition  takes  the  simple  degenerate  forms  34/3z  >  0  and  4  ~ 
respectively.  The  Green  functions  G®  and  (7"  corresponding  to  these  limiting 
cases  are  given  by  4xG®  “-l/r  -  l/r"*  and  4ir(?*  -  -1/r  +  l/r",  tdiere  r  -  l(C“x) 
+(n“y)^+(5~z)^ ]  «s  was  defined  previously  and  t* •  [  (5~x)^+(n-y)^+(j:+r)^] 

The  modified  Integral  equation,  corresponding  to  equation  (4),  for  the 
zero-frequency  potential,  4O  say,  takes  the  form^  ‘ 


4^(5)  -f0a)+L0(C;4°). 


where  the  potentials  f^(C)  and  1^(5 ;4®)  are  defined  as 


fO(g)  -  -/g  G0(5,x)40n(T)da(T), 


lO(T;40)  -  Is  G0„(t,T)[40(t)-40(T)]da(T). 


The  corresponding  oodlfled  integral  equation  for  the  infinite-frequency 
potential,  say,  takes  the  form 

[i-w(T)]r(T)  -  f(T)+L-(T;D.  (7) 

where  the  potentials  f*(C)  and  L**(C;^**)  are  defined  as 

f“(T)  “  “/s  Cr(C,x)^*„(3r)da^), 

-  Is  Gr„(t,t)[*-U)-*-(t)]da(T), 

and  w(^)  Is  the  waterplane  integral  defined  as 
2iiw(T)  -  -C  Ifa  [(e-x)^+(n-y)^+C^]“^/2dxdy, 

with  o  “  portion  of  the  mean  free-surface  plane  inside  the  body  (w  =  0  for  a 
fully-submerged  body)*  By  expressing  the  Integrand  of  the  foregoing  waterplane 
integral  in  the  form  3(x-C)[(y-Ti)^+C^]‘'^[  vx-5)^+(y-n)^+C^J”^^^/5x  »nd  using 
Green's  ^eorem,  we  may  express  the  waterplane  integral  in  the  form  of  the 
equivalent  waterline  integral 

2xw(c)  -  /c  (x-C)[(y-n)^+C^]"^[(x"€)^+(y-n)^+C^]“^^^dy, 

where  c  represents  the  intersection  curve  of  the  body  surface  S  with  the  plane 
z  “  0. 

FLOW  ABOUT  AN  ELLIPSOID  IN  AN  UNBOUNDED  FLUID 
The  practicality  of  solving  the  classical  integral  equation  (1)  by  using 
an  Iterative  solution  procedure  based  on  the  Neumann  series 

♦  (n+l)(T)/2  -  f(T)+L(T;^<">),  (8) 

with  n>0  and  =  0,  has  been  Investigated  numerically  for  the  simple  case  of 

flow  about  an  ellipsoid  with  major,  intermediate  and  minor  dimensions  equal  to 
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1,.  0.15  and  0.1.  respectively.  Figure  1  depicts  the  relative  errors  associated 
with  the  nt>»  approximations  the  added-mass  coefficients  for  all 

six  fundamental  motions  of  the  ellipsoid  (i  -  1  to  6).  This  figure  shows  that 
the  Neumann  series  (8)  converges  fairly  rapidly  for  all  motions  except  for  surge 
(translation  in  the  direction  of  the  major  axis).  This  numerical  result  is  in 

7 

agreement  %ri.th  the  theoretical  result  obtained  in  Noblesse  and  Trlantafyllou 
for  translatory  motions  of  ellipsoids.  In  particular,  it  is  proved  in  this 
study  that  the  Neumann  series  (8)  fails  to  converge  for  a  needle-like  spheroid 
in  longitudinal  translation.  Figure  1  shows  that  the  sixth  and  fourth  itera¬ 
tions  are  essentially  exact  for  sway  and  heave  (translations  in  the  directions 
of  the  intermediate  and  minor  axes),  respectively.  Convergence  is  even  faster 
for  roll  and  pitch,  for  which  three  iterations  are  sufficient.  Finally,  five 
iterations  are  necessary  for  yaw. 

The  velocity  potential  ^(C)  associated  with  the  modified  Integral  equation 
(4)  and  defined  as 

^»(T)  =  f^(T)/[f(T)-L'(T;f)] 

Is  proved  in  [?]  to  be  exact  for  translatory  motions  of  any  ellipsoid.  This 
theoretical  result  has  been  verified  numerically  for  the  ellipsoid  considered 
in  Figure  1.  The  potential  i|i(C)  has  also  been  verified  numerically  to  be  exact 
for  rotational  motions  (roll,  pitch  and  yaw).  The  accuracy  of  the  calculated 
added  mass  coefficients  for  all  six  fundamental  motions  is  within  0.4Z  by  using 
approximately  400  panels.  The  above-defined  potential  <|»(C)  corresponds  to  the 
first  approximation  in  the  sequence  of  iterative  approximations  4i(")(C) 

defined  by  the  recurrence  relation^ ^ 

♦  (n+l)(T)  -  ♦<")(T)f(?)/[^(">(T)-L'(T;^(">)], 

with  n>0  and  ^^^)(^)  =  f(C).  associated  with  the  modified  integral  equation 
(4).  The  practicality  of  using  this  recurrence  relation  for  solving  the 


modified  integral  equations  (6)  and  (7)  is  tested  numerically  in  the  following 
section  for  a  typical  hull  form. 

FLOW  ABOUT  A  SHIP  FORM  AT  ZERO  AND  INFINITE  FREQUENCY 

The  foregoing  recurrence  relation  for  potential  flow  in  an  unbounded  fluid 
takes  the  form 

^0(n+l)(O  .  ♦o<n)(T)fo(t)/[*0^”^<T)-Lo(T;tO^"^)]  (9) 

for  flow  in  the  zero-frequency  limit,  c^’-responding  to  the  modified  integral 
equation  (6).  The  corresponding  recurrence  relation  associated  with  the 
modified  Integral  equation  (7)  for  flow  in  the  infinite-frequency  limit  is 

The  recurrence  relations  (9)  and  (10),  with  n  >  0  and  “  fy,  ^ 

have  been  used  for  calculating  the  six  added-mass  coefficients  Aj^^  for  i  =  1  to 
6  in  the  zero-  and  infinite-frequency  limits  for  the  ship  form  depicted  in  Figure 
2,  with  beam/length  and  draft/length  ratios  equal  to  0,15  and  0.05,  respec¬ 
tively.  The  relative  errors  associated  with  the  n^^^  iterative  approximations 
A£i(n)  to  the  added  mass  coefficients  Ajj  for  all  six  fundamental  motions  of  the 
ship  (1  =  1  to  6)  are  depicted  in  Figures  3  and  A  for  flow  in  the  zero-  and 
infinite-frequency  limits,  respectively.  These  figures  show  that  six  to  eight 
Iterations  are  required. 

CONCLUSION 

The  numerical  results  presented  in  this  report  Indicate  that  the  recur¬ 
rence  relations  for  free-surface  flows^*^  corresponding  to  the  recurrence 


relations  (9)  and  (10)  for  flow  in  the  zero-  and  infinite-frequency  limits  may 
provide  an  efficient  solution  procedure  compared  to  that  associated  with  the 
classical  integral  equation  as  shown  in  Figure  1,  especially  for  cases  (hull 
form  equlped  with  a  bulb  or  a  sonar  dome,  low  Froude  number,  high  frequency)  for 
which  a  large  number  of  panels  is  required^ .  Further  numerical  study  of  these 
recurrence  relations  is  required  for  ascertaining  their  practical  usefulness. 
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Figure  1  -  Convergence  of  the  Neumann  Series  (8)  Associated  with 
the  Integral  Equation  (1)  for  Flow  About  An  Ellipsoid 


Figure  3  -  Convergence  of  the  Recurrence  Relation  (9)  Associated  with 
the  Integral  Equation  (6)  for  Zero-Frequency  Flow  About  a 

Typical  Ship  Form 
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THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS.  A  SEMIFORMAL  SERIES.  CONTAIN  INFORMATION  OF  A  PRELIM¬ 
INARY,  TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN¬ 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE-BY-CASE 
BASIS. 


